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Abstract. We consider eigenfunctions of the Laplace-Beltrami operator on special 
surfaces of revolution. For this separable system, the nodal domains of the (real) 
eigenfunctions form a checker-board pattern, and their number is proportional 
to the product of the angular and the "surface" quantum numbers. Arranging the 
wave functions by increasing values of the Laplace-Beltrami spectrum, we obtain the 
nodal sequence, whose statistical properties we study. In particular we investigate 
the distribution of the normalized counts ^ for sequences of eigenfunctions with 
K < n < K + Aif where K, AK e N. We show that the distribution approaches 
a limit as K,AK oo (the classical limit), and study the leading corrections in the 
semi-classical limit. With this information, we derive the central result of this work: 
the nodal sequence of a mirror-symmetric surface is sufficient to uniquely determine 
its shape (modulo scaling). 



1. Introduction 

Nodal domains of a real, continuous function are the maximally connected domains 
where the function does not change its sign. The nodal domains of eigenfunctions of the 
Laplacian on compact domains have been studied since Chladni first observed the nodal 
structure of vibration modes of thin plates (eigenfunctions of the bi-harmonic operator), 
in the early years of the 19**^ century. In the present manuscript we are not interested 
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in the geometric properties of nodal domains of Laplacian eigenf unctions, but rather in 
their count. 

Following Courant, we order the eigenfunctions so that the corresponding 
eigenvalues form a non- decreasing sequence. Denoting by Vn the number of nodal 
domains of the n-th eigenfunction, we form the normalized nodal sequence '■= ^i n & 
N. Courant 's theorem |16j guarantees that < 1? and we would like to study the 
distribution of the values of ^„ in the unit interval (0, 1]. In previous papers [H E] the 
distribution of the ^„ for various planar domains and 2-manifolds were studied, and 
it was concluded that the features of the distribution depend crucially on the type of 
classical dynamics it supports. If the classical dynamics on the manifold (geodesies) are 
integrable (and quantum mechanically separable - for such systems, actually, quantum 
separability is equivalent to classical separability t29j), the limit distribution exists, and 
displays certain features which are common to all such systems. On the other hand, if 
the classical dynamics is chaotic, the distribution of the normalized nodal sequence is 
well reproduced by using a random wave model for the eigenfunctions [5j. Bogomolny 
and Schmit computed the mean and the distribution by using ideas from percolation 
theory [1]. 

The first work with implications on the geometric content of the nodal sequence was 
of Smilansky and Sankaranarayanan [2], where it was shown that the aspect ratio of a 
rectangular domain on the plane (with Dirichlet boundary conditions) can be determined 
by counting its nodal domains. In [6] the nodal sequence for eigenfunctions of the 
Laplace-Beltrami operator for "simple" surfaces of revolution was discussed. A trace 
formula for the nodal count was derived, and was shown to depend explicitly both on 
some mean geometric properties of the surface, as well as the lengths of its geodesies. 
In spite of the formal similarity between the spectral and the nodal trace formulae, 
the geometrical information is included in different ways. Further studies [S, [TJ [32] 
have shown that isospectral domains have different nodal sequences, thus supporting 
the conjecture that the geometrical information is stored in the nodal and spectral 
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sequences in different ways. In the present work we go one step further and inquire 
whether one can deduce the shape (up to scahng) of a domain given the distribution of 
the normalized number of nodal domains, or, paraphrasing the classical spectral inversion 
question posed by M. Kac |19j, can one count the shape of a drum? It must be stressed 
that the only use of the spectral information is lexicographical - it is ordered as a 
non- decreasing sequence. Otherwise, there is no reference to the actual values of the 
eigenvalues. 

In the present manuscript we shall confine ourselves to the integrable (and 
separable) case, particularly to a special class of surfaces of revolution. We present 
new results which pertain to the distribution P(^, Ik) of the normalized nodal counts 
of eigenfunctions with indices n in the interval Ik = [K-, K + AK]. In [1] it was showed 
that there exists a limit distribution P{C) when the size of the index interval Ik becomes 
infinite (corresponding to the semi-classical limit, K oo). We provide the leading 
order term of the difference between P(^, Ik) and the limit distribution P{C)'- 



We show that the knowledge of the function P{^) and Pi{C) suffices for nodal domain 
inversion, provided that the surface is mirror-symmetric. In other words, given the 
normalized nodal sequence, we can deduce uniquely the profile function of the surface 
of revolution (provided it is smooth and symmetric). Numerical simulations were carried 
out for ellipsoids of revolution, which illustrate our theoretical findings. 

In what follows we shall use the classical notation of asymptotic analysis; the 
standard 'O, o' order notation, the symbol '~' standing for an asymptotic relation, 
the symbol 'x' denoting same order of magnitude, and the symbol denoting greater 
order of magnitude. 

2. Surfaces of revolution 




(1) 



We consider surfaces of revolution M. in R'^, which are generated by the complete 
rotation of the line y = f{x), x E I := [—1,1], about the x axis. We confine our 
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attention to a special subset of functions which satisfy the following requirements: 

i. f^ is analytic in /, and vanishes at ±1, where f{x) ~ a± {l^xy/^, with a± > 0. This 
requirement guarantees that M. is compact, has no boundary and is smooth even at the 
points where it is intersected by the axis of rotation. 

ii. The second derivative of / is strictly negative, so that f{x) has a single maximum at 
some X = Xmax-i where it reaches the value fmax- This requirement guarantees convexity 
of M. 

Surfaces which satisfy the requirements above will be referred to as simple surfaces 
of revolution, and are convex, mild deformations of ellipsoids of revolution. The induced 
Riemannian metric on Al is 

ds^ = (1 + f\xf)dx^ + f{xfde'' , (2) 

where the prime denotes differentiation with respect to x, and 6 e [0, 27r) is the azimuthal 
angle. 

In the proceeding subsections we shall review the properties of geodesies (classical 
mechanics) and the spectrum of the Laplace-Beltrami operator (quantum mechanics) 
on M. 

2.1. The geodesies 

The geodesies on M are the classical trajectories of free motion. They can be derived 
from the Euler-Lagrange variation principle with the Lagrangian 

L = 1 ((1 + nxf)x' + fixfe-") , (3) 

where a dot above denotes time derivative. The angular momentum along the axis 
of rotation f{xY9 is conserved, and we shall denote its value by m. The momentum 

conjugate to x is Px = (1 + /'(a;)^)x, and the conserved energy is 

2 

m 

It is convenient to introduce the action variable n, 



E^{1 + f{xf)x' + ^ . (4) 



n{E, m) i- j pAx = ^ £ ^/EfW^^^^^S^dx . (5) 
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Here, x± are the classical turning points, where Ef{xY — m'^ = 0, with x_ < x^ax < 2;+, 
which correspond to two meridians 7-1- (projections of caustics onto Ai) between which 
all geodesies with m 7^ wind around A4. Real classical trajectories exist only if 
E > {m/ fmaxY- The convexity of guarantees that the action variables {n,m) along 
with their conjugate angle variables constitute a global coordinate system on phase 
space [25] . 

The classical Hamiltonian H{n,m) in the action-angle representation is obtained 
by inverting ([5]) to express the energy in terms of n and m. H{n,m) is a homogeneous 
function of order 2, i.e. H{Xn,\m) = \^H{n,m), A > [13j. It suffices, therefore, to 
study the function n{m) := n(l,m), which defines a smooth line F in the (n, m)-plane 
(the projection of the unit energy shell on the action plane). The function n{m) is one 
of the main building blocks of the semi-classical theory which will be used throughout 
this work. We shall list some of its properties which will be used in the sequel: 

1. The refiection symmetry, n{—m) = n{m), follows from the definition (Q. Thus, we 
restrict our attention to m > when referring to F. 

2. n{\m\) is defined in the interval = (0, rrimax], where rUmax = fmax- In this interval, 
n{m) is analytic and decreases monotonically since 

dn(m) _ m f + 1 v^TTTW < ^ 



(7) 



dm IT J,_ ^f{xf - m2 f{x) 

3. n{m) assumes its maximum value at m = , 

n(0) = - j Jl + /'(x)2 dx = -, 

TT TT 

where C is the length of the rotating line. We show in Appendix A that n{m) is not 



analytic at m = 0, and in that vicinity 



n[m) ~ \m\ 

TT 



(8) 



At the other endpoint, n{m) vanishes. 



n{m) \l -{rrimax -m) , = l^nimaxf" {Xmax)\- (9) 



UJ 

4- The phase space volume is 



(2vr) 



e(E-H{n,m))dv = 2E n{m)dm = 2AE 

} 'J *i 



(10) 
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A is the area enclosed between the hne F and the n and m axes. It is related to the 
area oi M hy \\M\ \ = SnA. 

5. The computation of the higher derivatives of n(m) cannot proceed simply by taking 
the derivatives of (EI) - the resulting integrals diverge. To overcome this difficulty the 



integral defining n{m) needs regularization [2?j. This is done in Appendix A 



6. Some authors (e.g., [13j) prefer to use the Clairaut integral T instead of the angular 
momentum. They are related by 

X=^. (11) 

The twist condition is introduced in [13] to distinguish the class of simple surfaces of 
revolution, for which the dynamics are particularly simple. In the present notation, the 
twist condition is expressed by the requirement 

' 7^ , < |m| < m^ax (12) 

i.e. n{m) is either convex or concave on (0,mmax)- 

Throughout this paper we shall use the ellipsoid of revolution to illustrate 
graphically our findings. The ellipsoids are generated by /(x)^ = £^(1 — x^), with 
£ > being the eccentricity of the generating semi-ellipse. 

The action variable n(m) reduces to 

n{m) — 



1 - (I - s^)xlt^dt , (13) 

71 ./n 1 — 



where x± = ±a/1 — m^/e^ are the classical turning points. This integral can be 
evaluated in terms of elliptic functions giving, 



n[m) 



&E(1 - it) - (1 - e')m'K{l - -) - e'U{l -—,,!--) , (14) 



b b b - 

where 6 = + (1 — e'^)m? , and 

K(fc) := / [l-ksm^ey^l^dd, 
Jo 

E{k) := / (i-ksm^ey/^e 

Jo 

n(A;,/):= / {i~ksm^e)-\i-isin^e)-^'^de 

Jo 
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are the complete elliptic integrals of first, second and third kind respectively. 
Figure 1. shows the functions F := n(m) for a few ellipsoids of revolution 
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Figure 1. The curves F for ellipsoids of revolution with eccentricities e = 0.5 (prolate), 
e = 1 (a sphere) and e = 2 (oblate) are shown as dashed line, solid line and sparsely 
dashed line resp. 



2.2. The Laplace- Beltrami operator 

The Laplace-Beltrami operator on M. reads 



1 d f{x)d 1 ^^^^ 



f{x)(7{x) dx a(x) dx /(x)^ 86"^ 



where a{x) := a/1 + f'{xY. The domain of A are ^ G 14^2^ (/ x S^), required to be 
27r-periodic in 9. Under these conditions, the operator is self-adjoint, and its spectrum 
is discrete and non-negative. A is separable, and the eigenf unctions can be written in 
the product form '^{x, 6) = exp{im9) ipm{x), where m G Z. It is convenient to introduce 
a new variable t, through 

dt = ^dx , (16) 
fix) 



+ ( - E^,rrJ{x{t) f ) 



^n,m(t) = . (17) 
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which maps the interval [—1, 1] to MUjcxo}. For any m and eigenvalue En^m, the spectral 
equation for (fT5l) reduces to the Sturm-Liouville ODE 

"dt2 

The spectrum of the Laplace-Beltrami operator is doubly degenerate for all m 7^ 0, 
with En-m = En^m- The semi-classical spectrum is constructed by using the Einstein- 
Brillouin-Keller approximation [13], 

E'^% = H{n + ]^,m), ri G No, m G Z , (18) 

where H{n, m) is the classical Hamiltonian defined in terms of the action variables. 
The semi-classical approximation for the spectral sequence with m = assumes a very 
simple form. Since 

n + l = ^ [' Vl + f'ixydx = , (19) 

the semi-classical quantization condition reads: 



jpscl 
K,0 



7r(n + i) 



2 



^ J , n G No . (20) 

Because of the degeneracy of the spectrum, we have to choose a particular 
representation of the wave functions. We do this by associating cos{m6) with m > 
and sin(m6') for m < 0, i.e. 

{sm{m9) if m < 
(21) 
cos{m9) if m > . 

The nodal pattern of \1' is that of a checkerboard, typical for separable systems (as 
a matter of fact, the nodal pattern remains a checkerboard for any linear combination 
of the basis functions. It is only rotated around the symmetry axis of the surface). 
For m = 0, the number of nodal domains is z/„^o = n + 1, and for all other m, 
t'n.m = 2(n + l)|m|. In summary 

+ l)(2|m|+5„,o) • (22) 

To end this section we illustrate its content by an application to the simplest surface 
- the sphere, considered here as a surface of revolution with /(x)^ = 1 — x"^. 
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The action variable ^ can be computed explicitly 




m 



(23) 



Thus, H{n,m) = {n + |?«|)^, and the EBK quantization for the spectrum is -En,m = 
(/ + ~ /(/ + 1) where l:=n+\m\. The {21 + l)-fold de generacy follows immediately 
by counting the number of integer pairs {n,m) which satisfy I = n + \m\. 

Turning to the quantum description, the variable t defined in ( |T6l) can be explicitly 
computed, x = tanht. Writing E = l{l + 1), I G No, transforms ( ITTl) to 



which is equivalent to the Legendre equation. Finally, the number of nodal domains of 
spherical harmonics is known, and coincides with (122|) when the identification I = n+\m\ 
is made. 

3. Counting nodal domains 

We shall start this section by reviewing some of the general definitions and results 
obtained in pj, where the limit distribution of the normalized nodal counts was first 
studied. We shall then derive the next to leading term, and show that it provides further 
information on the geometry of M.. 

The nodal structure of the wave functions was reviewed in the preceding section, 
and an explicit expression for the dependence of the number of nodal domains z/„ ^ on 
the quantum numbers (n, m) is given in fl22|) . 

The object which is investigated in this work is the nodal sequence which is defined 
as follows: Arrange the spectrum as a non- decreasing sequence. This amounts to 
assigning to each pair of quantum numbers (n, m) a counting index A/'(n, m), which gives 
the number of eigenvalues of the Laplace-Beltrami operator (counted with multiplicity) 
which are strictly smaller than the eigenvalue -En,m; i-e. N'{n,m) := ij^{E G Spec(A) : 
E < En^m}- Obviously Af{n,m) = N{En^m), where N{E) is the spectral counting 
function. 




(24) 
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To account for the spectral degeneracy, we modify the definition above for |m| 7^ 0, 
so that ^/{n, \m\) = ^/{n, —\m\) + 1. The nodal sequence is the sequence of nodal 
counts ordered by Af: {i^m}m=i- ^^^^ convention, the systematic degeneracy of 
the spectrum is taken care of. In general, however, accidental degeneracies cannot be 
excluded. The ordering ambiguity may appear, e.g., when the degeneracy class involves 
states with different (non-negative) m values, such as e.g., for the sphere. In this case 
a possible way to remove this problem is to consider the sphere as a limiting case of an 
ellipsoid with (a positive) eccentricity approaching zero. The m degeneracy is removed 
for any ellipsoid, and the order of the eigenvalues is monotonic in |m| for arbitrary small 
eccentricity. Similar constructions can be used for other accidental degeneracies. 

Courant's theorem [1^ ensures that, for any ordering of the eigenfunctions in the 
degeneracy classes 

J^Af(n,m) < A/'(n, m) . (25) 

It is natural therefore, to define the normalized nodal sequence 

:= ^ , < ^AT < 1 . (26) 

We study the distribution of the values of the normalized nodal sequence for a finite 
index set, Af e {K,...,K + AK}, 

^(e, Ik) := ^ Yl '^)) - ' (27) 

(n,m,)6NoxZ 

and xik the characteristic function of the interval Ik, 



1 if K <x<K + AK 
XI Ax) = { (28) 
otherwise . 

Since -P(^, Ik) is not a function but a distribution over the interval (0, 1), one must 
take care in its manipulations. Most of the limits and estimates are considered in the 
weak sense (e.g. [E]). In some sections though, related distributions are viewed as 
functions (or to be more precise, the functions whose values these distributions take on 
some subinterval of (0, 1)). In other cases they will be manipulated as functions after 
an appropriate regularization. 
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In PP, it was assumed for convenience that Ik grows linearly in K, AK = gK, g 
being a positive constant (such that gK G N) . The existence of the limiting distribution 
P(^) of P{^, Ik) in the K oo limit (the classical limit) was proved, and its universal 
features were presented. The existence can be proven by showing 



AK 



PiO^m , (29) 



uniformly in g, for any smooth and compactly supported test function (f on (0, 1). Here, 
we shall repeat the derivation in more detail, and also compute the difference between 
P{^,Ik) and P{i) to leading order in (again, for the linear case AK = gK - the 
superlinear case AK ^ K, follows trivially). 

We rewrite fl271) as a sum of an isotropic component, to which only isotropic states 
(with m = 0) contribute, and an anisotropic component (with m ^ 0), 

P((Jk) =f„..o(f,-f/{) + i'„#oK,/if) (30) 

p^Mi, - ^ E X.. mn, »))*(«- . 

Since we are interested in the semi-classical limit, we are allowed to make the following 
approximate steps, which incur errors of order higher than 0{^^). 
i. The spectrum is approximated by -E„,m ~ E^'^j^ = H{n + |,m), as m? + n? oo, 
which introduces a relative error bounded by 0{^—) [12]. In particular, as was shown 
in (EOD, = . 

a. To the same order, A/'(n, m) can be replaced by the first term in the Weyl series |13j . 

Af{n, m) = 2AH{n + ^, m) (l + 0(^)) , (31) 

2 V 7714 / 



-'n,m 



where A was defined in (fTOj) . 

Introducing these approximations in fl5U]) . we find that Pm=o{C,, Ik) is 0{^) (in 
the weak sense) and therefore we defer its computation to a later stage. The sums 
over [n, m) in fl30l) are computed using the Poisson summation formula, decomposing 
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-Pm^o(^) ^k) in a smooth and an oscillatory part, 

Pn.M^, Ik) = Ik) + Q{^, Ik) ■= P{^, Ik) + Q n,m{^,Ik) , (32) 

(Ar,M)eZ2\{0} 

where the Fourier coefficients are 
Qn,m{^, Ik) ~ 

/r *" /.I <''~'X,. i2AHin + 1 - 5i|^)dn . (33) 

and obviously P(^, J^) = Qo,o(^, ^k)- 

The leading term in the sum above is the smooth term P{C,, Ik) which we calculate 
first. The oscillatory terms are of lower order in and will be computed in a separate 
subsection. 

Proceeding with the smooth part, we shift the integration variable n i-^ n + | and 

write the result as 

2 / \ / (n + -)m\ 

P(^M~^f^i,nf ,,,^(2^H(,,™))i(^-l^)dn . (34) 

We change the integration variables {n,m) i— > {S,s), where S = H{n,m), and s is 

defined through the relations 

, , dH(n,m) dH(n,m) 

at = UJndn + CiJmdm ; UJn = , UJm ~ 



dl — III r^ 1 

n o m 

ds = YT~2 • (35) 

Note that with this definition the Jacobian is unity, and dndm = d£^ds. Thus, ( 134|) is 

reduced to 

K(l+g) 

where G(x) is the Heavyside step function. The pair of functions {n{s),m{s)} constitute 
a parametric representation of the line F, along which H{n{s),ni{s)^ = 1. This allows 
the scaling by which appears in fl36l) . 

The expression above can be further reduced by the following observation. On the 
line F we have Undn + uJmdm = 0, which induces a symplectic structure (where n and 

m play the roles of canonically conjugate variables and s is the "time") 

dm(s) _ _ dn{s) _ 
as ds 
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Another important identity follows from the fact that H{n, m) is homogeneous of order 
2 in {n,m). We can write n = n{H,m) = -^n{l, -^), from which we deduce that on 
the line {H{n,m) = 1}, 

— -LOn - —, 7 -, r , [6^] 

as n[m) — mn'[m) 



dn(m) 

. y(l+g) 



where n'{m) := . Thus, fl36|) takes the form 



ZlA 57? V / 

A similar transformation can also be applied in the computation of the oscillatory 
integrals Qn,m{^^^k)- 



3.1. The limit distribution 

The limit distribution is the leading term in the expansion of the integral fl36|) in powers 
0[^^). Taking the hmit K — oo, the £^-integral can be directly performed, resulting 
with 

P(^) = \l iU-^^\is. (40) 



This expression allows a simple geometrical interpretation: the product n{s)m{s) is the 
area of the rectangle whose vertices are the points (m(s),n(s)) on F, its projections 
(m(s),0) and (0,n(s)) on the two axes and the origin (See Figure 2). P{Cj is the 
probability distribution of the areas of these rectangles (scaled by A., and therefore 
smaller than 1). The areas of rectangles which are based on points which are either near 
the m or the n axes approach . Since F is either convex or concave, there exists a 
unique point where ^ reaches the maximal value of scaled areas, ^max- Thus, -P(0 = 

for i > irnax- 

Changing integration variables using (!38|) . we get 

P(^)= / 5 U ^]\n{m)-mn'{m)\dm . (41) 



2A Jo V A 
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Figure 2. The geometric interpretation of P(f) (Computed for an ellipsoid of 
revolution with eccentricity e = 0.5). 



Integrating (jH} we get, 



IE 



n{m)—mn'{m) 
n{m)+mn' {m) 



f- mn(m) 



(42) 



The sum is over the real values of m > which satisfy ^ = in the vicinity of 

imax, typically two solutions coalesce, leading to a square root singularity of -P(0 that 
point. The vanishing of -P(0 iii the interval [^max,!], and the square root singularity 
at ^max are the universal features which characterize the nodal domain distributions for 
separable systems (in 2-d) in general, and simple surfaces of revolution, in particular. 
Using dm), it is easy to check normalization, P(^)d^ = 1. 

The form fH2|) for the distribution can be further simplified since for simple 
surfaces of revolution the twist condition f|T2|) is satisfied, and there are only two values 
of m which solve ^ = VUl^nl^ Denote them by m^(^) and m^(^) (m_(^) < m+(^)). 
They coalesce at ^ = ^max- The value of m where the sole maximum of mn{m) occurs 
is denoted by mo = m±{^max)- The two functions m+(^) and m_(^), together, provide 
a parametric representation of the curve F, since n±[^) = where < < ^max- 

This parametrization will be used often in the subsequent discussion. 
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To leading order, m±{^) ~ "^0(1 ± a/TCI)! where ( 
deduce that, in the left neighborhood of ^rnax 

max 



and so from 

<^m,ax 
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we 



(43) 




Figure 3. P{£,) for ellipsoids of revolution with eccentricities e = 0.5, 1 and 2 (sparsely 
dashed line, solid line and dashed line resp.) 

One can re-arrange (I42p to obtain another expression for the limiting distribution. 

For ^ < ^maxi 

This expression is quite revealing, because it can be inverted to provide the function 
based on information derived only from the nodal sequence. We shall show below 
that the next to leading expression in the expansion of P{^, Ik) provides another relation 
between m_(^) and m^{C,). Solving the two equation we obtain a complete parametric 
representation of the curve F (or n{m)). We shall also show that F defines uniquely 
the function f{x), when f{x) is symmetric about x = 0. This will prove our claim that 
the nodal sequence for symmetric surfaces of revolution can be inverted and provide the 
"shape of the drum" ! 

The behavior of P{C,) near ^ = can be easily extracted using (jH]). Since ^ —>■ 
implies either m — > or n{m) —>■ along F, we may use the linear approximations 
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for n(m) as given by ([HD and Q), respectively. The equation ^"^^ = ^ reduces to a 
quadratic equation, and its solutions define the two branches and as 



271 \ V " £2 



Substituting in (jHl) we get 
where, 



The relation (146!) shows that limg^o+-P(0 = ^1 independently of the surface under 
consideration - another universal feature to be added to the aforementioned ones. 
Moreover, it shows that one can extract the dimensionless geometric parameters ri_ and 
?7+ from PiX) in the neighborhood of = 0. They are directly related to the properties 
of the line which generates M. through its length, maximum distance from the axis of 
revolution and its curvature at the maximum. Note, however, that the nodal sequence 
is composed of integers, and in contrast to the spectral sequence it is invariant under 
isotropic scalings on M.. Therefore, only dimensionless quantities can be extracted from 
it. Here, all the lengths are expressed in units of \fA. (The parameter r/o in fH7j) is 
another dimensionless parameter which we define here even though it will appear only 
later) . 

The limit distributions for three ellipsoids of revolution are shown in Figure 3. 
Computing ^max as a function of the eccentricity reveals that imax{^^ is a monotonically 
decreasing function which varies between ^maa'(O) ~ .550 and ^maa(oo) ~ .464 as shown 
in Figure 4. Thus, one can deduce the eccentricity from the nodal sequence just by 
determining the support of -P(0- 
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0.55 



to 



0.5 



0.45 








4 



Figure 4. £,max as a monotonically decreasing function of the eccentricity e for 
ellipsoids of revolution. The dashed line marks e = 1 - the sphere - for which £,max = -5 
as in ([481) . 

As another application of (jS]), and as an illustration, we derive -P(0 for the sphere. 
From (1A.12P we get that for the sphere n(m) = 1 — m which immediately gives 



The same result can also be obtained directly. The spectrum consists of the values 
En,m = n{n+l) which are independent of m (|m| < n) and are (2n + l)-fold degenerate. 
The eigenfunctions are the spherical harmonics, and for the sake of counting their nodal 
domains, we consider them in their separable basis. The number of zeros of P^{cos6) is 
n when m = 0, and n—\m\+2 otherwise. The number of nodal domains are respectively 
n + 1 and 2|m|(?2 — |m| + 1). The counting function for the 2n + 1 degenerate states 
with a given n satisfies 



n 



(49) 



The ordering within the set is arbitrary, but in the limit of large n it is immaterial, since 
to leading order we can write M{n^m) = n{n + 1)(1 + 0{^)). Hence, the contribution 
of the n-fold degenerate set to -P(0 is 



n(n + 1) 



+ 0i-) 
n 
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j\(i-2^l{l-^^)y^^ 



1 



18 

(50) 



3.2. The next to leading terms 

Several terms contribute 0{^^) corrections. In the following we shall address them in 
detail. 

The leading correction to P{^,Ik)- Starting again from ( l39ll . we see that the m- 
integration results in a sum of two terms, coming from the two branches of solutions 
^ — ^ '|^(™-) _|_ ^j^^ _|_ (9(1) which constitute the two-point support of the S function 
in the integral. They differ by amount Sm±{^) of 0{-^) from the values m±{C,) which 
were introduced for the computation (jUj) of the limit distribution, 
1 m 



5m±(0 = — ^ — — + 

2y 8 ran [m) + n[rn) t 



computed at m = m±{^) . (51) 



As long as (m + 5m) ± remain inside the integration range, one can proceed with the 
computation towards the semi-classical asymptotic expansion 



K, 



P 



Of course, P is the limit distribution (H2l) . and 

pm = 



9 



where p±{C) ^-re given by 



n(m) 



1 + m 



n'(m) 
n{m) 



1 — m 



n'(m) 
n{m) 



1 + m 



n'(m) 
n{m) 



-1 + 



p_(0+P+(0 



2m ( m " j-"^} 

n(m,) 



(52) 



(53) 



m 



n'(m) ^2 
n{m) ' 



computed at m = m±{^) .(54) 



These expressions can be further simplified. To do so, we differentiate rn±{^)n{m±{^)) 
A^, with respect to ^ and get 



n 



m 



n 



n 



1 



dm 



1 

m 



n" = n 



(55) 



d5 '-""d? '-Vd? 

computed at m = m±{^). By substituting the above in the defining expressions for 



P±(0) obtain after some straightforward calculations 



I ( dm± d2m±\ 



dm-t 



(56) 
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Hence, 



A(0 




2 



9 




This is an explicit expression which provides the leading correction in terms of the 
difference — between the two branches of the parametric representation 

of r. Together with it forms the basis for the inversion procedure which will be 
discussed in detail in the next section. 

The conditions for the validity of the above approximation are not satisfied if ^ 
is in the O(^) neighborhood of either or ^max- Near ^max (15T!) diverges, while 
near ^ = 0, (m + 5m)± may lie outside of the integration range To get 

the correct expressions for -P(^, Ik) in the vicinity of the extreme values of ^, we use 
several variations of the same trick: within the problematic domains of integration, we 
approximate n{m) as a linear function of m. The argument of the 5 functions become 
quadratic functions of m. The support of the 5 functions can be evaluated explicitly, 
and the m-integrations can be performed exactly. The remaining £^-integrations turn 
out to be straightforward, so that explicit expressions for P(^, Ik) in the vicinity of 
^ = and ^ = ^max are obtained. 

Behavior of P{^, Ik) in the neighborhood of ^ = ^max- We recall that ^max is defined 
as the maximum value of HL-^Uli^ which occurs at mo, where n{mQ) + mo n'{mQ) = 0. 
Thus, in the neighborhood of mo we can approximate n{m) ~ no — — "^o) where 

no = n^mo) = '^^^ (the non-existence of other critical points is guaranteed by the 
twist codition). With this approximation, the argument of the 5 function in ( 139|) is 
quadratic in m. The integrations over m and £ have to be carried out with attention to 
the requirement that the support of the 6 remains within the integration range. After 
some lengthy but straightforward manipulations one gets. 




(58) 



and 



(4no)2 



C 



max 



(59) 



'^max 



2A 



max 
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Vmax is another dimensionless parameter which characterizes P{^, Ik) near the borders 
of its support. Performing the integral we get, 

^(e, Ik) ~ riT^ ( Fiy^^ - Fiy^, C) j (60) 




=arctanh 



—7^ arctan 



K\+y 



^/\C\^/y''+'2y-<: 
-(+y 



VCy/y^+2y-( 



C < 



1 



Kri„ 



Kr]„ 



9)KVn 



C 



max <-,.,,, 



for C < 



for C > . 



This expression includes the corrections to the limit distribution at its most noticeable 
feature, namely its singularity at C,max- For finite K, the square-root singularity is 
replaced by a continuous function which reaches beyond C,max, shifts the maximum from 

^max to ^maxi^I^^d) 

e 



and extends the support of P(^, Ik) up to 



+ 2 



KTJmax 



(l+g)Krimax 

As ^ oo, the expression converges to the limit. The application of 



the above theory for an ellipsoid of revolution for two values of K are shown in Figures 
5. and 6. 
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Figure 5. Local behavior of P(^, Ix) &t ^ — £,max for an ellipsoid of revolution with 
eccentricity e = 0.5 {K = 24000 and g = 1). 
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Figure 6. Local behavior of P(^, Ik) at ^ = S^max for an ellipsoid of revolution with 
eccentricity e = 0.5 [K = 48000 and g = 0.5). 



Behavior of P{^, Ik) in the neighborhood of ^ = 0. In the neighborhood of = we 
have contributions from the support of the delta function in (139|) from the neighborhood 
of m = - provided that (m + 6m) ^ > - and m = rrimax - provided that 

^ > -^fmn(m) + • Putting all contributions together, we have 

\ / TfLrn.n.T. 



H^Jk) 




1 1 



9 

1 



:i+9)-^ 



1 

+ 2< 



1 



9 

1 



[1 + 9)-^ 



for 

for 

for 
for 
for 
for 



^(1+9) 



K <^ 

<e< 



K{l+g) 



VP 



Ka+9) 



K{l+9) 



(61) 



The most important feature in fl^T]) is that it shifts the support of Ik) away from 
^ = to min < . . k Note that most of the expressions which make up 



P(^, Ik) in the neighborhood of ^ = are confined to a ^ interval of size 0{^^). Beyond 
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this interval (-^ <^ ^ < ^max) Ik) takes the form 

+ , (62) 




which coincides with the small ^ expression of the limit distribution fj46l) to leading 
order in ^. To obtain the dominant behavior of P{^, Ik) near ^ = 0, one should add 
Pm=o{^, Ik) which we compute now. 

The contribution of the m = term in ^3^). Following the same steps as above, the 
leading approximation to Pm=o{C,, Ik) reads, 

. _L J- (2Amn + i. 0))) * - dn. (63) 

This being already a correction term, we are allowed to neglect the semi-classical 
correction ^ to n in the argument of H{n + |,0). With H{n,0) = (^)^ from (120]) . 
we get 

/ K(l + g~ 

1 /•vV^:4- 1 £2 



2A 



p...„«.w~^y M?--^jd„. (64) 




e gK \l Kil+g) V K ^ ^gg^ 

otherwise 
where //_ was defined in (l47j) . Both the shift of the support away from the origin ^ = 
and its size, decrease as In its support, Pm=o{C,, Ik) is bounded between the values 
i and 1 + ^- However, although the contribution of this term to the probability density is 
0(1), its effect on the probability is 0{^), or in other words, Pm=o{^, Ik) = 0{^) in the 
weak sense. In this vicinity, Pm=o{^, Ik) depends only on a single geometric parameter, 
ri_. Note that a priori, one would have expected the lower limit of the support of 
P{C,,Ik) to be 0{^), and not 0{^^). This prediction stems from the definition of 
the normalized nodal sequence (126|1 : for A/" > 1 we have the lower bound ^j\f > jf, so 
P(^, Ik) = for ^ < jj^^, and not for < 0{-^) as observed. 

The oscillatory terms Q{Ik) contribute terms of order in the C, = vicinity. 



They originate from a Gibbs phenomenon and they will be discussed in [Appendix B 
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Figures 6. and 7. compare the results of numerical simulations with the expressions 
derived above in the vicinity of = 0, and the theory includes also the oscillatory 
corrections (to be precise, with the integrated density n(^,//^) := P(^', /;^)d^'). 




0.01 0.02 0.03 0.04 0.05 0.06 



Figure 7. Local behavior of n(^,//f) at ^ = for an ellipsoid of revolution with 
eccentricity e = 0.5 [K = 24000, g = 1). 




Figure 8. Local behavior of n(^,//f) at ^ = for an ellipsoid of revolution with 
eccentricity e = 0.5 {K = 48000, g = 0.5). 
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3.3. The oscillatory contributions Q{^,Ik) 

Weak estimate of Q{Ik)- In this section we estimate the oscillatory part Q{^,Ik), 
defined in ( 132|) . We show that it is of smaller order in than the leading order 
correction of the corresponding smooth part. This justifies the preceding analysis where 
we considered only the smooth part, which gives the only contribution of 0{-^). The 
numerical results also give evidence of this self-averaging process in the semi-classical 
limit. 

To be more precise, we shall show that 

Q(e,/x)^(Ode = o(-i) , (66) 

for any test function (p. We shall follow a series of natural regularizations which are 
justified in studying this weak limit. We shall also discuss the local behavior of Q{^, Ik) 
in the neighborhood of ,^ = 0, in order to compare with the numerical results. 

We have Q(e, Ik) = E (Af,M)^o QN,Mi^j ^k), where the Fourier components ([33]) are 
approximated as 

Qn,m{^, Ik) ~ 



1 r{'i-+9)K/2A rmmax mr)(m\ n- 

' d£ 5{^- !^^!^) nim) - mn'{m) e2-v^(*^™+^«(™))dm . (67) 



9K Jk/2A Jo ^ 

We now turn to a smoothing of this distribution by adding a small imaginary part 

to the argument of the delta function, say e > 0. This amounts to replacing the S 

function by the Lorentzian 

1 

And so, 

^ l'{l+g)K/2A ^2niy/£{Mm+Nn{m)) 

Qn,aM^ ^k) 7- / -TTTTc — 7: \^ ~ mn'{m) dm . (69) 

T^gK JK/2A ^ U(m),m) 



5^(0 = 23fJ / e2"'«"e-2"""da; = — . — . (6^ 

' — _L 



We must notice here a qualitative difference between the semiclassical theory of the 
spectral density and the nodal domain distribution. The index set Ik = [K-, K + /S.K] , 
corresponds to a spectral interval [Ek, Ek+i^k]- In the later K — » 00, the 

eigenvalues are distributed in an ever growing interval, while the normalized nodal 
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sequence is distributed in (0,1], becoming arbitrarily dense. So, in contrast to the 
regularization of the spectral density, in general we do not expect the limits £ 0+ 
and i^T — > oo to commute for the nodal domain distribution. We shall first consider the 
semi-classical limit. 

Once again, by the linear approximation, we have n(m) = n±+n'^{m—m±)+0{^), 
and mn{m)/A = ^{n± + m±n'_^){m — m±) + 0{^) =: ^±(m — m±) + 0{^). 

By taking the whole real line as the m-integration range and shifting m ^ x — 
m — 777.0- {(^ — we have 

}K/2A foc JlmVeRx 

<K/2A J 



rdx 



E 



1 n„-m^< r{l+9)K/2A 2^.5^_2^^^|« 



^ gKna + m^n'^ Jk/2A 
2 n.-mX /■V(i+5)^/2'4 



—r^ / e ' kdk , 

gK + m^n'^ J,/WJU 



(70) 



with 



R n' 



\ ^ J \^'^ J 
Now, instead of a sharp uniform window (taking into account only those states 

whose index lie in the interval Ik) , we consider a Gaussian regularization, 

hK{k) = exp ( - (A; - f,)yv') (71) 



■^-Y^y Ta ^ ~ (-\/l + g — ^)\ r^i and extend the integration over the 



with - 2 \l 2A 

whole real line (of course, this regularization does not affect the asymptotic behavior 

we study), 

2 

<T=± 



= E 



n„ - m„n„ 



yfiivgK Ua + rrian'^ 

gK Ua + rrian'^ 



-k'^/v^+(2-n:iS-2Tre\^\+2n/v^)k, 

e ' €a ' ' ' kdk — 



iriS — ne 



R 



/i 



X 
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X exp f - n'^v^S'^ + vr^e^t;^^ - 2n^iev'^ t- ^ + '^l^-'^iS - 2ith6 

We proceed towards a crude, yet sufficient for our purpose, estimate for Q{^, Ik 

2v^ 



R 

Sct 



26 

.(72) 



gK 



Tin 



Ha + m„n' 



Ti\S\ + Tie: 



R 

So- 



+ 4 'X 

V 



{N,M)y^O cr=± 

-j^ - 2nixe - 

By taking into account that "^^^^2 e~^^'^"^''^^^^"^ x ^ and X]7ez2 |7*N 
t —>■ oo, we have 



X exp ^ - n^v'^S'^ + Tx'^e^v^^ - 2nfxe ^ j 



(73) 

t(a-7)2+v^/3-7 ^ 



|Q.(e,/i.)l<0(-^), (74) 

uniformly in e (before the hmit e 0+ is taken, all three terms in the amplitude of the 
summed quantity contribute to this same order; x -^)- 

Local estimate of Q{C,, Ik) the neighborhood of C, =0. From the numerical 
investigations, a clear oscillatory behavior of the distribution is observed in the 
neighborhood of ^ = 0, which implies the importance of the oscillatory part in that 
region. This is due to a Gibbs phenomenon, and is not in contradiction to the fact that 
the oscillatory part Q{Ik) is of less order in the weak sense than the smooth part P{Ik) 
- this is a local contribution. 

We find the manifestation in the dominant terms X^atQa^.o and ^jvf Qo,j\/, which 
agrees with the numerical results. The explicit form of these contributions are presented 



m 



Appendix B 



4. The geometric information stored in the nodal sequence 

The derivation of the probability density -P(^, Ik) in the preceding section is based on 
functions which are obtained directly from the dynamical relations embedded in the 
Hamiltonian and its dependence on the action variables. These, in turn are computed 
from the profile curve y = f{x), which defines Ai. Here, we would like to investigate 
the possibility of inverting this relationship, and ask what can be said about the surface 
once the nodal sequence is known. 
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A few parameters can be easily extracted from the probability density P{^,Ix)- 
Taking the limit ^ oo we get the limit distribution P{^). Its support provides 
the parameter C,max- In the vicinity of ^ = 0, P{C,) depends symmetrically on the two 
parameters and 77+. This provides a useful relationship, which can be combined 
with the information from P{^, Ik) near = 0, for finite K, which depends on all four 
(independent) geometric parameters, rj_,rj^,rjQ and rj^ax- These parameters are related 
to the geometry of the surface (H71) . 

As was commented in the previous section, the relations (jH]) and (1571) . can be 
solved to obtain the two branches m^{^) and m_|_(^). Consequently, the function n(m) 
is determined, since n{m±{^)) = j^^^^- This will be the first step towards the inversion 
of the nodal counting data which will be explained in the next subsection. 

4.I. Inversion of the Nodal Sequence 

In this section we discuss nodal domain inversion in detail. The first result is: 
with the distribution P{Ik) Cis given data, the "scale-invariant" action n{m)/\fA is 
determined uniquely. For simple surfaces of revolution (separable systems in general 
with the Hamiltonian satisfying the homogeneity condition), n{m) incorporates all the 
information about the dynamics, and as we shall see in some cases, the actual geometry. 
Our starting point is the semi-classical asymptotic expansion (j52l) . 




(75) 




(76) 



Obviously, P and Pi are determined from the given data. 



lim P{Ik) = P , and 




(77) 



From relation , 




(78) 
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the ratio rn^{^)/m_{^) is determined, given P. Accompanying the above, we have the 
relation (|5711 . 



(79) 



which can be inverted to give 

uiO ^(™.K)-™-«))^-^^^/^(/^dr)dr(SO) 

To determine m(^) uniquely, two integration constants should be provided, and they are 
given in terms of the initial conditions for u{C,) and ^ at the point ^ = or ^ = ^max- 
They are expressed in terms of the geometric parameters {r^., 77+, ?7o, ?7maa}, which are 
extracted from the form of -P(0 near ^ = and ^ = C,max- Thus the ratio and the 
difference between m^(^) are m_{^) given. Together with n±{^) = ^ ^^(-g-, they give the 
parametric representation of F. Note that we require no information from the statistical 
properties of the nodal counts of the isotropic quantum states - i.e. Pm=o{lK)- 

The significance of this result becomes apparent in the next section, where we 
confine ourselves to mirror-symmetric surfaces. 

4-2. Mirror- symmetric f{x) is uniquely determined by the n{m) 

Following the preceding section, we prove that if the generating curve / is mirror- 
symmetric, (i.e. f{—x) = f{x)), the action variable n{m) determines the / uniquely. 
Since / is an even function, we may write 



2 



TT 



x+ 



n[m) = — I \/f[xY — m'^— — — — ax 







/(^) 



where /(x^) = m, or x+ = / ^{m). Changing the variable x ^ u = f{x)/mr, 



Suppose that /, g are different curves of this class which give the same action 
variable, i.e. 



^/rrimaxifyu'^ - m"^ + f O f ^immax{f)u)'^ 

(/)«)! 
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Jm/mma.{g) 19 ° 9 \m^ax{9)u)\ 

while f 9- 

Firstly, note that rrimaxif ) = ^max{,9)- This is because nimax is the sole real root of 
n{m), thus, since both sides of the above equality are proportional to n{m) they must 
have a common root, denoted simply be rrimax, since n(m) is analytic on {0,mjnax]- 
Thus, since the integration limits are identical, the integrands must be equal. 

It suffices to show that = ^^^Xg-Hvlf '' ^ ^ [^^rrimax], implies 

f = 9- From the above we have /' o f~^{y) = 9' o 9^^{y)- The problem has been 
reduced to showing that the nonlinear operator Af = f o f^^ acting on our function 
class, possesses an inverse, i.e. Af = Ag ^ f = g. 

Consider the inhomogeneous "functional" equation f of~^ = h, for some h in some 
other appropriate function class. We shall show that for given h, this determines, along 
with some initial condition, a unique /, formally / = A~^h. We have f'{f^^{y)) = h{y), 
or f'{x) = h{f{x)), since y = f{x). Thus, we have reduced this to a ffist order ordinary 
differential equation, y' = h{y). Accompanied by the initial condition /(O) = mmax, 
this becomes an initial value problem on [0, 1] with a unique solution {h is smooth on 
[0,1)). Thus, A-^ exists. 

In our problem the initial condition is provided by knowledge of the root of n{m), 
so we have reached the conclusion that the two integrals cannot equal ii f ^ g. 

5. Summary and conclusions 

The unique inversion of the nodal sequence which was demonstrated above for symmetric 
and "simple" surfaces of revolution paves the way to a sequence of problems which 
should now be addressed. Other families of separable manifolds are known, amongst 
which the Liouville surfaces [30] and the axially symmetric Zoll surfaces [3T] are of prime 
importance. We believe that the general approach taken in the present paper could be 
applied to handle these case, however, modifications should be applied to take care of 
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special problems which are intrinsic to these problems, and this remains for a further 
study. Another class of integrable systems consists of independent particle models which 
are commonly used in Atomic and Nuclear physics. The study of such systems extends 
the research of nodal domains to systems with arbitrary dimensions. 

The next systems in complexity are systems which are classically integrable but 
are not separable quantum mechanically. The simplest examples consist of e.g. the 
Dirichlet Laplacians in the equilateral or the isosceles right triangles. Even though the 
spectrum can be expressed precisely in terms of the "quantum numbers", counting of 
nodal domains is difficult, and the study of the nodal sequences in such cases might call 
for other approaches then the one pursued here. 

Do nodal sequences in other systems store geometric information? Is there a way to 
extract this information to determine the geometry? These are yet open problems, and 
the only hint for an affirmative answer comes from preliminary numerical simulations 
which indicate that "nodal" trace formulae exist for "quantum graphs" [32] and for 
"chaotic billiards" [33]. No rigorous treatment exists so far. 
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Appendix A. The regularization of the action integral 

The action variable ([5]) is defined in terms of an integral, whose form is not convenient 
for further computations, such as e.g., the evaluation of its higher derivatives. This can 
be done by regularizing the integral in a way which will be explained here. Rather than 
introducing fractional derivatives as was done in e.g., [23], we compute the integrals 
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explicitly. It is convenient to introduce the notation q{x) = /(x)^ and /i = m^. We 
start with 

, , 1 r+ V4g(a;) + q'{xy , , , , ,^ ^, 

^lAi) = - / Vli^) - A* ?r-r-\ ' where q{x±) = ^ . (A.l) 

TT A_ 2g(x) 

The function g(a;) is analytic in / and has a single maximum at Xmax- We separate the 
integration interval in flA.ll) to two consecutive intervals [x^,Xmax] and [xmaxja^+J and 
write accordingly 

n{fi) = n^{fi) + n+{fi) . (A. 2) 

In each of the intervals [—1, Xmax) and {xmax, 1], lix) is a monotonic function. Therefore 
it can be inverted in each of the intervals in terms of the corresponding functions x±{q), 
which are analytic in the interval [0,qmax)- Thus, 



1 



Qinax 



n±{yi) = - Vq - /iw±('?)dg (A. 3) 

71 



i.Lfi^^iMy+i. (A.4) 



with 



2qV \ dq 

The expressions in the square brackets above are analytic in the domain of q where x±{q) 
are analytic, that is, in [0, qmax)- They can be Taylor expanded with a convergence radius 
Qmax SO that 

^ r=0 

Substituting this expressions in ( 1A.3I) . performing the integrals and defining = 
r+ + r ~ , we finally get 



71 



{qmax - /i)^ - arccos 



oo 

+ - Vr, , (A.6) 



r=0 



where, 

= g Vq^dq = {qmax - /^) ^ 2^ (^^J _ ^) ^ 3 • (A-7) 

Using ([7]) we find 
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Thus in the vicinity of /i = we get 

C 

n{fi) v^. (A.9) 

TT 

The behavior of n(/i) near the other extreme end of the interval - /i = Qmax - cannot be 
deduced in the same way, because x{q) is not defined at this point. However, staring 
directly from flA.l|) we can obtain the behavior of n{fi) in this domain. For this purpose 



we write 



q{x) ~ Qmax -7;^{X- Xmaxf , = \q" {,Xmax)\ , X Xmax ■ (A.IO) 



2 

To leading order in {qmax — Z^), ( 1A.1I) reduces to 



^ /O / \ = /o ~ \ -Umax - \m\)- (A.ll) 

7rV2a;gmax V ?max - q ^/2ujqraax V 

Finally, we return to the example of the sphere. With f{xY = 1 — x"^ the integrals 

can be performed exactly, 

n(m) = 1 — \m\ . (A. 12) 

This is consistent with the local expressions presented in flA.9|A.lip . 



Appendix B. The oscillatory part near the origin 

We recover the oscillatory behavior of the distribution near = from the leading order 
sums 

Q{iJK) ~ QNflii^K) + QoM^^Ik) . (B.l) 

We begin with ([321) • As ^ 0+, the contribution in the above integral will come from 
the neighborhouds of m = and m = rrimax (the delta function of the integrand is 
supported on {m + 5m}±). These will be treated separately, denoted by ]^{^, Ik) 
respectively, so that 

Qn,m{^-> ^k) = Qn^mI^^ ^k) + QN,Mi^^ ^k) ■ (B.2) 

In what follows, we define x := ^/rj^ when referring to the Q~ terms, and x := C/r]+ for 
the terms. 
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Given that m_ ~ -^x and 5m_ ~ —\tix^ by changing variables to the appropriate 
dimensionless wavenumber k :— jCVE, we have 



<5o,m(C, Ik) ~ 

Similarly, we carry out the calculation of the integral Ik)- Here, the condition 

that (m + 5m) lies in the integration range reads > ( mn{m) + r]m ) = 'mmaxV: 



TTlma 



1 e" 



^T^M^/^{l+x/2) r^Jm{l+g)K 



VoK 



VW2 



r]ogK 1-x 
where the appropriate wavenumber is A; := rUmaxV^- 
Following the above calculations, 

Qn,o{CJk) ~ 

1 p,mNx/2 /■7r-^??-(l+g)K 

TT^rj-gK 1-x J^,J^ \ xJ V 2/ 



dA; . (B.4) 



and 



Q1,A^,Ik)^ 

1-x //^ \ xJ \ 2 J 



TjogK ±-x jvw^ 
By performing the integrations, we have 





Mez* 



1 



1 



1 



for 



nrj^gK (1 — a;)a;^ 



fo^ a/2(IT^>^ 



T:r]ogK (1 - x)(2 - x)^ 



where 



Ik) for ^ 2r,o(i+5)K < ^ ^ 



(B.5) 



(B.6) 



fii^^lK) for ^^j^-i-^<e<y^ 



X, 
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/^-(e, Ik) a{x^/{l + g)rJ^K - |) - a{x^f^ - |) + 



+ xV^/9(x^/^- |) + 9)r]-K^{x^{\ + ^)7?_i^- |) , (B.7) 



and 



/k) ^= - 1 + (1 + f ) ) - - + (1 + 1)^1)+ 



:=a((2-a;)V(l + ^)^oi^ + (l + |)y|)-«((2-a;)V^+(l + ^)y|' 

+2V;^/3((2 - x) (1 + \)'^- 

- 2V(1 + ^)r7oi^/?((2 - x) V(l + 5)^oi^ + (1 + 1)^1) ■ (B 



Similarly, 



5] QAr,o(e,^K) 



1 







for 



7rrj_gK (1 - a;)(2 - x)^ 



^1 ^ir) for 



1 



(1+9)V-K 



^ (e, Ik) for < C 



a; 



for 



2(1+9)%/^ 



27rriogK (1 — 



where 



2 X 
9i (C, ^x) := a(- - 1 + -) - a I 



(B.9) 



.((2-x)v^+|) + 

^^-(e, ^i.) «(| + (2 - x)y/{l + g)r^^K) - a(| + (2 - x) + 
+ 2 ((2 2^/(1 + ^)^-^/3 ((2 -x)V(l + ^)r;-i^ + ^) , (B.IO) 
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and 

jJK, W := a{.^{l + gf-^K- (1 + |)) - a{x^j^ - (1 + |)) + 

Above, we have made use of the periodic functions a and /3, which we define via 
their Fourier series, 

, , 1 e'^'"* 2 cos Trnt , ^, 

and 



neZ* nePJ 



e'^'"* v^sinTrni ,^ 
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